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We study primordial gravitational waves from inflation in the Randall- Sundrum braneworld 
model. The effect of small change of the Hubble parameter during inflation is investigated us- 
ing a toy model given by connecting two de Sitter branes. We analyze the power spectrum of 
the final zero-mode gravitons, which is generated from the vacuum fluctuations of both the initial 
Kaluza-Klein modes and the zero mode. The amplitude of fluctuations is confirmed to agree with the 
four-dimensional one at low energies, whereas it is enhanced due to the normalization factor of the 
zero mode at high energies. We show that the five-dimensional spectrum can be well approximated 
by applying a simple mapping to the four-dimensional fluctuation amplitude. 
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I. INTRODUCTION 



Our four-dimensional universe might be embedded as a three-brane in a higher dimensional spacetime with all the 
standard model particles confined to the brane and with gravity allowed to propagate in the extra dimensions. This 
so-called braneworld scenario opened various possibilities. A possible solution to the hierarchy problem in particle 
physics was presented by introducing large extra dimensions 1] or a new type of compactification . Also interesting 
is a new possibility proposed in Ref. 3] that the four-dimensional gravity is recovered effectively on the brane despite 
the infinite extension of the extra dimension 0, Q . 

The braneworld scenarios have also had a large impact on cosmology (for a review of cosmological aspects of the 
braneworld scenarios, see, e.g., JjJ). Homogeneous and isotropic cosmological models have been built Q, and various 
types of inflation models proposed d, 0, [H], 0> E3 ■ A lot of effort has been put forth in searching for new and 
characteristic features of the braneworld cosmology. For example, cosmological perturbations and the related physics 
of the early universe have attracted large attention in the expectation that the braneworld inflation might leave their 
characteristic prints on the primordial spectrum of perturbations. 

Alth oug h the cosmological perturbations have been discussed in a number of publications (see, e.g., |tl ITTl IT^ . Il3l 
li^lllillall^llll ll^ and references therein), the presence of the extra dimension does not allow detailed predictions 
about the cosmological consequences. As the simplest case, the primordial spectra of density perturbations and of 
gravitational (tensor) perturbations were investigated neglecting the nontrivial evolution of perturbations in the bulk 
in a model of slow-roll inflation driven by a inflaton field confined to the brane 0, fl3l Il4j . (We must mention 
that gravitational waves have been considered in the context of braneworld models other than the Randall-Sundrum 
type as well; see |2fJ,|21j.) Although complicated, the effect on the scalar type perturbations due to the evolution of 
perturbations in the bulk has been discussed under some assumptions 01 ■ As for the gravitational wave perturbations, 
the authors of Ref. considered a simplified inflation model in which de Sitter stage of inflation is instantaneously 
connected to Minkowski space. In this model it is possible to solve the perturbation equations including the bulk to 
some extent. They focused on perturbations with comoving scale exceeding the Hubble scale at the end of inflation. 

Also in this paper we study gravitational waves from an inflating brane. If the Hubble parameter on the brane is 
constant, the power spectrum becomes scale-invariant |13|. However, the Hubble parameter usually changes even dur- 
ing inflation. The change of the Hubble parameter, i.e., the nontrivial motion of the brane in the five-dimensional bulk, 
"disturbs" the graviton wave function. As a result, zero-mode gravitons, which correspond to the four-dimensional 
gravitational waves, are created from vacuum fluctuations in the Kaluza-Klein modes as well as in the zero mode. It 
is also possible that gravitons initially in the zero mode escape into the extra dimension as 'dark radiation'. Therefore 
we expect that the nontrivial motion of the brane may leave characteristic features of the braneworld inflation. If so, 
it is interesting to search for a signature of the extra dimension left on the primordial spectrum. However, there is a 
technical difficulty. When the Hubble parameter is time dependent, the bulk equations are no longer separable. Then 
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we have to solve a complicated partial differential equation. To cope with this difficulty, we consider a simple model 
in which two de Sitter branes are joined at a certain time; namely, we assume that the Hubble parameter changes 
discontinuously. In this model we can calculate the power spectrum almost analytically. This is a milder version of 
the transition described in [l8| . 

This paper is organized as follows. In the next section we describe the setup of our five-dimensional model, and 
explain the formalism introduced in Ref. to solve the mode functions for gravitational wave perturbations. Using 
this formalism, we explicitly evaluate the Bogoliubov coefficients in Section 3. In Section 4 we translate the results 
for the Bogoliubov coefficients into the power spectrum of gravitational waves, and its properties are discussed. We 
show that the power spectrum for our five-dimensional model can be reproduced with good accuracy from that for 
the corresponding four-dimensional model by applying a simple mapping. Section 5 is devoted to conclusion. 



II. DE SITTER BRANE AND GRAVITATIONAL WAVE PERTURBATIONS 

A. Background metric, gravitational wave perturbations, and mode functions 

Let us start with the simple case in which the background is given by a pure de Sitter brane in AdSs bulk 
spacetime. We solve the five-dimensional Einstein equations for gravitational wave perturbations. For this purpose, 
it is convenient to use a coordinate system in which the position of the brane becomes a constant coordinate surface. 
In such a coordinate system the background metric is written as 



ds 2 = 



(sinh£) 2 



(2.1) 



where I is the bulk curvature radius, and the de Sitter brane is placed at £ = const = £g. Note that here r\ is supposed 
to be negative. On the brane, the scale factor is given by a(rj) — l/(—r)H) and the Hubble parameter becomes 

H = e^smh^B- (2.2) 

Note that under the coordinate transformations 

t = i] cosh £ — rjo cosh^B, 

z = — 7ysinh£, (2-3) 

with a constant 770 , the metric (|2.1|) becomes the AdSs metric in the Poincare coordinates. 
The metric with gravitational wave perturbations is written as 



dS = (sTnW ^ ~ ^ + {SlJ + W)****! + * >■ 12 ' ! 



We decompose the transverse-traceless tensor hj^ into the spatial Fourier modes as 



W(V,*,0 - ■ ^372 J d "P oul-PU-'-,,,. (2.5) 

where is the polarization tensor, and the summation over different polarization was suppressed. M5 represents the 
five-dimensional Planck mass, and it is related to the four-dimensional Planck mass Mp\ by £Mf = Mpj. The factor 
(Afs) 3 / 2 is chosen so that the effective action for <fi corresponds to the action for the canonically normalized scalar 
field. Then, the Einstein equations for the gravitational wave perturbations reduce to the Klein-Gordon equation for 
a massless scalar field in AdSs, 

□0=(2>„- D f ) 4> = 0, (2.6) 

where the derivative operators are defined by 

1-1 2 U o 1 2 2 

V v=V -rr^ - 2rj—+p 77 , 

^ = (sinh0 3 ^(sinhO- 3 ^. (2.7) 



3 



We assume Z 2 symmetry across the brane. Assuming that anisotropic stress is zero on the brane, Israel's junction 
condition gives the boundary condition for the perturbations as 



(2.8) 



Since the equation is separable, the mode functions are found in the form of </> K (?7,£) = V , k(??)Xk(C)i where ip K (v) 
and Xk(0 satisfy 



(2.9) 
(2.10) 



T>£ + K + ~) x «(0 = 0, 



respectively, and the separation constant n (> 0) is related to the Kaluza-Klein mass m as m 2 = (k 2 + 9/A)H 2 . For 



-9/4, we have one discrete mode, which is called the zero mode. The zero-mode wave function is given by 



(2.11) 



which is independent of £. The factor C{H) is to be determined by the normalization condition, (<j)^ ■ <^>o ) — -Flj 



with respect to the Klein-Gordon inner product 

(F ■ G) := -2i f 

Then we have 



£ 3 d£ 



(sinh^) 3 77 ; 



(Fd v G* - G*d v F) 



(2.12) 



C 2 (H) 



2(sinh^) 2 / d£ 



(sinh£) 3 



y/l + £ 2 H 2 + i 2 H 2 In 



IE 



,1 + s/l + PH 2 , 

This normalization factor is the same that was introduced in, for example, Refs. [l3l Il9l ] , and behaves like 



(2.13) 



C 2 (H) 



1, (at low energies £H <C 1), 

3 

~£H, (at high energies IE > 1) 



(2.14) 



The continuous spectrum called Kaluza-Klein (KK) modes starts with k — 0. Notice that the mode labeled by 
k = does not correspond to the zero mode @ . Writing the positive and the negative frequency modes, respectively, 

as </>l + ^ = ipK Xk an d <pk = (0k + ^)* = iprT^Xni we impose the conditions 



il 3 



(sinhO 3 



X* k 'Xk = S(k-k'), 



(2.15) 
(2.16) 



so as to satisfy the normalization condition (</>L ± ' ) • <j>^ ) = T<5( K — The solutions of Eq. I|2.9J) are given in terms 
of the Hankel functions by 



The spatial mode function Xk (£,) is given in terms of the associated Legendre functions by |llj 



(2.17) 



Xk = dismht;) 2 Prr /a+te (coshO - C 2 QZi /2+iK (coshO 



(2.18) 
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where from Eqs. (|2.8() and (|2.16(l the constants C\ and C 2 are 



Ci 



Co 



r(— i/c) 



r(5/2 + m) 

-Pri /2+tK (cosh^) 
<9li/ 2+iK (coshes)' 



r(5/2-i«) 



ttCi 



r(i/t-3/2) 



r(i + ««) 



n-l/2 



(2.19) 



As will be seen, we need to evaluate the value of the wave function at the location of the brane Xk(£b), and in some 
special cases Xk(£b) reduces to a rather simple form. For sinh£s <C 1 and Ksinh^s < 1, we have 



k tanh ttk / k 2 + 1/4 
2 \ k 2 + 9/4 



(sinh^) 2 



while, for sinh^s > 1 or Ksinh^s 3> 1, we have 



1 



K 



Xk (Zb) « ^(sinhC S ) 3/2 - _ . 

V7T ^K 2 +9/4 



(2.20) 



(2.21) 



For the derivation of these two expressions, see Ref. [l8[ . 



B. Model with a jump in the Hubble parameter 

We consider a model in which the Hubble parameter changes during inflation. As we have explained, in the case 
of constant Hubble parameter the brane can be placed at a constant coordinate plane. When the Hubble parameter 
varies, we need to consider a moving brane in the same coordinates. For simplicity, we consider the situation in which 
the Hubble parameter changes discontinuously at 77 = rj from Hi to 

H 2 = H X - SH. (2.22) 

Here 5H/H\ is assumed to be small. For later convenience, we define a small quantity e# by 



tHi y/i + (w 2 y - w 2 y/i + (my 

ch = 



IH 2 

2 / xu\ 2 /m\ 3 



I SH , 2 + 3(^0' (SHY + W\ (2 . 23) 



y/l + {my #1 2(l + ((iH 1 ) 2 Y' 2 \Hi) \H X 

To describe the motion of the de Sitter brane after transition, it is natural to introduce a new coordinate system 
(77, £) defined by 

t = fj cosh £ — 770 cosh £b , 

z = -7?smh|. (2.24) 

Then, the brane expanding with Hubble parameter H 2 is placed at £ = £b by choosing two constants £s and 770 so as 
to satisfy H 2 = i~ x sinh£s and 770 sinh^s = ijo sinh^s- The trajectory of the brane is shown in Fig. ^ Apparently, 
mode functions in this coordinate system take the same form as those in the previous section, but the arguments 
(£, rj) and the Hubble parameter Hi are replaced by (£,fj) and H 2 . We refer to these second set of modes as (f>o and 
<f> K . The relation between (77, £) and (fj,Q is 



fj = -\JV 2 + 2e ff 77o?7coshe + e 2 H r] 2 , 
tanh£ = (T/cosh^ + e#7i n ) _1 7isinh£. (2.25) 

As explained above, the variation of the Hubble parameter is assumed to be small. For a technical reason, we 
further impose a weak restriction that the wavelength of the perturbations concerned is larger than SH/H 2 . These 
conditions are summarized as follows: 

-Jf<h (2-26) 
SH 

p|r? |-^«l. (2.27) 
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FIG. 1: Trajectory of the brane (thick solid line) in static coordinates. Dashed (respectively dotted) lines represent surfaces of 
£ = constant (respectively £ = constant). 



C. Method to calculate Bogoliubov coefficients 



We consider the graviton wave function tp that becomes the zero mode 4> ' at the infinite future fj — 0. We write 
the wave function tp as 



<f = 0o ) + Sl P> 



(2.28) 



where the second term Sip arises because (f> ' does not satisfy the boundary condition for t < 0. Writing down ■* 

and Sip at the infinite past, r\ = — oo, as a linear combination of 4>q ^ and <f>^ , we can read the Bogoliubov coefficients. 

(-) 



At r\ — > — oo the first term 0jj is expanded as 



6^ — ► E ( ^4' + VbM^ ; , 

ri — > — oo * — ' V ' 
M=0.k 



(2.29) 



where the summation is taken over the zero mode and the KK modes. The coefficients Uq k and Vo« are written by 
using the inner product as 



Uom = lim 

Tj — > — OO 



Vqm 



lim 



V • A?) , 



77 — >•— 00 

Evaluation of the inner product at 77 = — oo leads to Vqo = Vqk = while 



u. 



oo 



H 2 C(H 2 ) ie H pT,a cosh^i 

#lC(-ffl) 



1 - i£Hp?7o (C 2 (i?i) - coshes) - -(e H p?7o) 2 (sinh^ B ) 



(2.30) 

(2.31) 
(2.32) 



These expressions are approximately correct as long as p\r)o\SH/Hi <C 1 is satisfied. The derivation of these equations 
is explained in Appendix iBl 

The second term Sip in Eq. H2.28[l is obtained as follows. Since both tp and cf>Q ' satisfy the Klein-Gordon equation 
in the bulk l|2.6l) . 5<p also obeys the same equation. The boundary condition for Sip is derived from 12. 8|) as d^Sip = 



-d, 



Wo 



Therefore, the solution Sip is given by 



no 



Sip = 2 I drfG ad v(ri,tW,t B ) %4 \rj',0 



(2.33) 
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with the aid of the advanced Green's function that satisfies 

(23, - 23 S )G adv (77, = S(rj - rf)5(£ - £')• 

The explicit form of the Green's function is 01 

^ (sinh^') 3 ?)' L 
Taking the limit 77 — > —00, we can expand Sip in terms of in-vacuum mode functions, 

Sip — 



[ u vm4>m ] 



Af=0,K 



where the coefficients are given by 



-2if 
vqm = 2U 3 



dr, 



(-)/ 



(sinh £,b) 3 V 



4^ fa 6*) 'faO 



(-)/ 



(sinh^ B ) 3 ?7 



(2.34) 
(2.35) 

(2.36) 

(2.37) 
(2.38) 



To evaluate these coefficients, we need the source term d%< 



written in terms of the coordinates (??,£), which 



is 



= -r z/2 sinh £ B G \-^ ipeHmV sinh ^e !p Vf 2 + 2e «*i rash ^+ e H''o. 



(2.39) 



From Eqs. I|2.29|) and (|2.36|l . we finally obtain the Bogoliubov coefficients relating the initial zero mode or KK 
modes to the final zero mode, 



<P 



M 



(2.40) 



where 



OtQM — Uqm+Uom, 

Pom = vqm- 



(2.41) 



From these coefficients we can evaluate the number and the power spectrum of the generated gravitons. 



III. EVALUATION OF BOGOLIUBOV COEFFICIENTS 



Now let us evaluate the expressions for the Bogoliubov coefficients obtained in the preceding section. We concentrate 
on the two limiting cases: the low energy regime {£H\ <C 1) and high energy regime (IH\ 3> 1). We first evaluate the 
coefficients aoo and /?00j which relate the initial zero mode to the final zero mode. We keep the terms up to second 
order in eu (or equivalently in 8H/H\). 

Substituting Eq. (|2~35|l into Eq. f2~3%|) . we have 

/3oo = C(Hi)C(H 2 & H r, I" dr, (\ - -U e M^Vif^ HV0VCOshiB+eWQ ) 
Hi J-oo \r) W) ) 

Because there is a factor e# in front of the integral, we can neglect the correction of 0(e 2 H ) in the integrand. Then 
we can carry out the integration to obtain 

/3oo « C(Hi)C(Hz)^-e H -^— e' 2tmo ^ motHCOsh ^ B . (3.2) 
Hi 2pr) 

Similarly, we get 

«00 « t/00 + C( J ff 1 )C( J ff 2 )^e H (\ + -?—) e -iPVoe H co S h( B ^ (3 3) 

Hi \ 2pr, ) 
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where Uqo is given by Eq. I|2.31[l . 

At low energies (£Hi C 1), the Bogoliubov coefficients ctoo and Poo become 



I i 5H 



2 OTo 



e -ipr)aSH/H 2 (3 4) 



/3 00 pa _J_?E e -2m,a-ipria5H/H 2 _ /35s 

2p?7o ffi 

It is worth noting that these expressions are correct up to the second order in 5H/Hi. This result agrees with the 
result of the four-dimensional calculation oj( 4D ) and P^ 4D K 
At high energies {£H\ ^ 1), the coefficients are 



1 x 3 ( i SH\ i 3 fSH\ ipr] 5H (pr] ) 2 fSH\ 



2\2pr]o Hi ) %\H X ) 2 H x 2 \Hi J 



e -ipVo ijf -ipva I ( I7) (-35) 



2 V2 P 77o #1 

Here we stress that the last two terms in the square brackets of Eq. (|3.t)|) . both coming from Uqq, are enhanced at 
PM > 1. 

Next, we calculate the Bogoliubov coefficients ao K and /3q k , which relate the initial KK modes to the final zero 
mode, up to the leading first order in 8H/H\. Although we will calculate the power spectrum up to second order in 
SH/Hi in the succeeding section, the expressions up to first order are sufficient for cyq k and fto K , in contrast to the 

case for a o and Poo- Again, substituting d^(p ~\r],^B) into Eq. Q2.38J1 . we have 

(3 0K = JTp 1 3/2 C(H 2 ) X UZ B ) S ' m ^ B e H7]Q P ^jr)(,) e W?+^^^. (3.8) 

(smh^s) 2 J-00 T 



Setting eh in the integrand to zero, the coefficient reduces to 



p\vo\ 



(3o,^ J ^CmxU^-^eHPVoe-^ 2 dxx- 3 / 2 Hg\x)e ix , (3.9) 



where we have introduced the integration variable x = —pr). Similarly we have 

OT-f rplvol 
^( J ff 2 )x:(C S )^ I ™oe-^ /2 y oo dxx-^ 2 H^(x)e- ix , (3.10) 

where U 0k is given by Eq. I|2.32|l and is C(e//)j2^]. 

Now let us discuss the dependence of a^ K and /3q k on £Hi and SH / H\ in the limiting cases, £H^ <C 1 and IHi ^> 1. 
At low energies, we find, using Eq. Q2.2(jp . 

|/?o.| 2 , |a 0K | 2 oc (IHJ 2 (j^j , (tHt « 1) , (3.11) 

where we have omitted the dependence on k and These coefficients are suppressed by the factors of £Hi and 

SH/Hi. Recall that aoo and Poo agree with the standard four-dimensional result at low energies. Thus, because of 
the suppression of ao K and (3q k at low energies, the four-dimensional result is recovered only by the contribution from 
the initial zero mode. 

At high energies, we obtain from Eq. I|2.21|l 

8H X 2 



\Po,\ 2 , |«o K | 2 oc , (Iff! » 1) , (3.12) 

where we have again omitted the dependence on k and p|f7o|- In contrast to the result in the low energy regime i|3.11|l . 
this high energy behavior is not associated with any suppression factor. 

Integrating |/?ok| 2 over the KK continuum, we obtain the total number of zero- mode gravitons created from the 
initial KK vacuum. It can be shown that the coefficients behave as Po K ~ ~ a 0K ~ (phol) 1 ^ 2 at p\rio\ <C 1, and we 
have 

/3 OK + a5«~(3(pho|) 3/2 . (3.13) 
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Thus the number of gravitons created is proportional to p outside the horizon and is evaluated as 



f 

Jo 



\f3 0K \ 2 d K : 



SHY 



Hi J 



(m « i), 



5H 



0.3 x p|i7o| { — ) , (£Efi»l) 



On the other hand, making use of the asymptotic form of the Hankel function H^\x) 
x —> co, we can evaluate the integral in Eq. I|3.9|) in the p\r]o \ — > co limit as 



i(x-(2v+l)n/A) 



dx- 



1 



P\Vo\ 



(p\Vo\ -> oo), 



(3.14) 



for 



(3.15) 



where we have carried out the integration by parts and kept the most dominant term. This shows that the creation 
of gravitons is suppressed well inside the horizon. Since the assumption of the instantaneous transition tends to 
overestimate particle production at large p , the number of particles created from the initial zero mode and KK 
modes is expected to be more suppressed inside the horizon than Eq. (|3.15(l if we consider a realistic situation in which 

the Hubble parameter changes smoothly. Note that uq k oc p\vo\ J^ V °^ x~ 2 dx is constant for p\r/o\ — > co. Therefore, 
the ao K coefficient is dominated by Uq k at large p, which behaves like |ao K | 2 ~ |£/o K | 2 K (PVo) 2 ■ 

An example of numerical calculation is shown in Fig. [2J The figure shows that the spectrum has a peak around 
the Hubble scale and then decreases inside the horizon, and we confirmed that the behavior of the number density 
outside the horizon is well described by Eq. 1)3 .14|) . 



log 
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FIG. 2: Spectra of number density of zero-mode gravitons created from vacuum fluctuations in the initial KK modes. Integration 
over k is performed numerically. Solid line represents the low energy case with £Hi — 0.1, while dashed line shows the high 
energy case. The latter should be understood as the limiting case £Hi — * oo since at high energies the number does not depend 
on IHv, see Eq. I5TT21 . 



IV. POWER SPECTRUM OF GENERATED GRAVITATIONAL WAVES 

So far, we have discussed the Bogoliubov coefficients to see the number of created gravitons. However, our main 
interest is in the power spectrum of gravitational waves because the meaning of "particle" is obscure at the super 
Hubble scale. 

A. Pure de Sitter brane 

Gravitational waves generated from pure de Sitter inflation on a brane have the scale invariant spectrum 

2C 2 (H) ( H 
'pi 



p -liTUl' (41) 
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which is defined by the expectation value of the squared amplitude of the vacuum fluctuation, 87rp 3 |</>o| 2 /(27rM5) 3 , 
evaluated at a late time; see Eq. (|2.11[l . Since C 2 ~ 1 at £H <C 1, this power spectrum agrees at low energies with 
the standard four-dimensional spectrum 

p » = 4© 2 ' (4 ' 2) 

At high energies, however, the power spectrum i|4.I[l is enhanced due to the factor C(H), and is much greater than the 
four-dimensional counterpart. This amplification effect was found in Ref. |l3j |. These results say that the difference 
between Eq. (|4.1|) and Eq. (|4.2[1 is absorbed by the transformation 

H i > HC{H). (4.3) 



B. Model with variation of Hubble parameter 



Now let us turn to the case in which the Hubble parameter is not constant. Time variation of the Hubble parameter 
during inflation brings a small modification to the spectrum, and the resulting spectrum depends on the wavelength. 
Here we consider the amplitude of vacuum fluctuation of the zero mode in the final state. It will be a relevant 
observable for the observers on the brane at a late epoch because the KK mode fluctuations at super Hubble scale 
rapidly decay in the expanding universe due to its four-dimensional effective mass. Since the behavior of the zero- 
mode wave function at the infinite future 77 — s- is known from Eq. i|2.11[) . we find that the vacuum fluctuation for 
zero mode at a late epoch is given by 



lim 



aoo</>o +) _ Poo4>o 



PoA ) 



C 2 (H 2 



Hi 



2p 3 



\aoo + Pool + / dn\ao K + Po K 



(4.4) 



Multiplying this by (2/M|)(p 3 /27r 2 ) and recalling the relation £M^ = M™, we obtain the power spectrum 

V 5D (P) = KS°(P)+-P^( P ), (4.5) 



with 



2C 2 {H 2 ) 


(H 2 


M 2 


\2w 


2C 2 (H 2 ) 


(H 2 


M 2 


\2w 



|«00 + Pi 



I In 



dK,\a 0K + Pq k \ 



(4.6) 



The power spectrum in the four-dimensional theory, computed in the same way, is given by 



M 2 



,(4D) 



■P 



,(4D) 



(4.7) 



The appearance of the coefficient a in the power spectrum may look unusual. This is due to our setup in which the 
final state of the universe is still inflating. In such a case, the fluctuations that have left the Hubble horizon never 
reenter it. Outside the horizon, the number of particle created, does not correspond to the power spectrum. 

There are two apparent differences between V5D and Vad ; the normalization factor C (H) and the contribution from 
In the low energy regime, however, the two spectra agree with each other: 



the KK modes Vf£ 



7>5D~7>4D («i«l). 



(4.8) 



This is because, as is seen from the discussion about the Bogoliubov coefficients in the preceding section, the zero-mode 
contribution 'Pfg' is exactly the same as Vati (up to the normalization factor), and the Kaluza-Klein contribution 
"P,^ is suppressed by the factor (£Hi) 2 . On the other hand, when £ti\ is large, the amplitude of gravitational waves 
deviates from the four-dimensional one owing to the amplification of the factor C{H). 

We have observed for the pure de Sitter inflation that the correspondence between the five-dimensional power 
spectrum and the four-dimensional one is realized by the map ()4.3|) . It is interesting to investigate whether the 
correspondence can be generalized to the present case. It seems natural to give the transformation in this case by 



h h-> hC(h), 



(4.9) 
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where 



h( P ) = h 2 



V (4D) 



=(4D) 



(4.10) 



namely, the rescaled power spectrum V rcs (p) is defined as 



2C 2 (h) ( h 



2tt 



(4.11) 



We will see that this transformation works well and mostly absorbs the difference between the five- and four- dimen- 
sional amplitudes. 

We examine the differences between the five-dimensional spectrum and the rescaled four-dimensional spectrum by 
expanding them with respect to 5H/ Hi as 



V 5D {p) 



v. 



(0) 



,(2) 



5D 



5D 



5D 



SHY 



H J 



O 



-jKK 

5D ' 



SHY 
~H ) 



(4.12) 
(4.13) 



Here the quantity associated with the superscript (n) represents the collection of the terms of 0((SH/H) n ). On 
the right hand side of Eq. I|4.12|) . all the terms except for the last one come from the initial zero mode. The direct 
expansion shows that the leading terms in Eqs. I|4.12(l and (|4.13|l exactly agree with each other up to the first order 
in SH/Hx, 



?(0) 

5D 



5D 



?(0) 
res 



v {1) 



2C 2 {H 1 



Hi 
2tt 



sin(2p?7o) 



PVo 



C 2 (iJi) 5H 
VI + (£H 1 ) 2 ~H~i 



(4.14) 



These terms contain only the contribution from the initial zero mode. 

The contribution from the initial KK modes is of order (SH/Hi) 2 because both ao K and /3q k are 0{8H/ H\). 
Thus, to examine whether the agreement of the spectrum continues to hold even after including the KK modes, we 

(2) (2) 

investigate the second order part of the spectrum. The second order terms and "Pros are given by 



2C 2 {H 1 ) fH, 



M 2 



2-k 

2f„„ \ r<2 



(PVo) 



C^jH^ + C 2 ^) 



-cos(pryo) 



C 2 (^) 



sin 2 ( Wo ) C 2 ^) sin(2pr) Q ) / 2 + 3(*ffi)' 



(PVo) 2 y/l + JJEcJ 
4C 2 (iJi) 3 + 4(£ffi) 



2pr] 



6C 2 (gi) 



Vi + (^i) 2 i + (^i) 



c 2 0ff! 



8_H_ 



2C 2 (H 1 ) (H x 



M 2 X 



2tt 



2cos(p77o) 
5(^) 2 



sin 2 (2pr to ) / 

(2 Pm ) 2 [ 1 + (M,) 2 V 1 + (^i) 5 
AC 2 (H^ 3 + 4(Wi 



s/i + (m) 2 

sm 2 ( CTo ) | sin(gj»ft) ( 2 + 3(lffi) s 

(pm) 2 pm 

4C 2 (ii 1 ) \ 



(4.15) 



(^i) 2 



4C 2 (ffx) 

VT+WhTY 2 



v/i + (^i) 2 i + (^i) 2 



C 2 (ffx 



Vi + (^i) 2 V^i 



(4.16) 



(2") (2) 

From these equations, we notice that Vrd and 7-5 D do not agree with each other. We see that the difference is 
enhanced in particular at p\r)o\ ^> 1. However, as we have mentioned earlier, the KK mode contribution also gives 
a correction of the same order, and in fact we show that approximate agreement is recovered even in this order by 
adding the KK mode contribution. 

(2) (2) 

First we observe the power spectrum at p\rjo\ <C 1. Expanding and 'P res with respect to p|?7o|, we have 



Vjlip) 
75(0) 







(4.17) 
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On the other hand, Eq. (|3.13(l leads to 



1 5D 



(P) 



J>(0) 



(p\vo\Y 



SHY 



(4.18) 



Therefore the difference is small outside the horizon as 



-p( 2 ) i T)KK 
' 5D r 5D 



T,(2) 
/ res 



p(o) 



(4.19) 



although the cancellation between P^(p) and Pi£l(p) does not happen. By a similar argument, at p\r]o\ < 1, we have 



,(2) 



,(2) 
5D 



•pKK 
' 5D 



(2) 



SHY 



(pM<i). 



(4.20) 



The situation is more interesting when we consider the spectrum inside the Hubble horizon. There is a term 
proportional to (pt]o) 2 in P^, which is dominant at p\rjo\ 3> 1, while there is no corresponding ter m in l^rd • Hence, 
the difference between and Vrel is 



? (2) 

5D 



Pi 



(2) 



-(PVo) 2 2- 



c- 2 (iJ 1 ) + c 2 (H 1 )\ c 2 (i?o 



\A + (^i) 2 / v/i + ^i) 2 



iff 



(4.21) 



Our approximation is valid for p\r)a\8H/Hi < 1 l|2.27|l . Hence, within the region of validity, this difference can be as 
large as . We also note that the terms proportional to (p\r]o\) 2 come from a o, while the contribution from /3 00 is 
suppressed at p\t)q\ ^> 1. On the other hand, the contribution from initial KK modes, V^^(p), is dominated by a.Q K 
at p\t)q\ 3> 1: P$q oc J dn\ao K + (3q k \ 2 ~ / d,K\ao K \ 2 - This means that, although the creation of zero-mode gravitons 
from the initial KK modes is negligible inside the horizon, a part of the amplitude of the final zero mode comes from 
the initial KK modes losing their KK momenta. Since the coefficient cto K is proportional to p\r]o\ at large p\r]o\, P^d 
behaves as 



Pf£(p) 
-p(o) 



iPVo)' 2 



SHY 



(4.22) 



(2) 

This KK mode contribution cancels P^ . The cancellation can be proved by looking at the property of the Bogoliubov 
coefficients 



I "oo I 



dK(\a 0K \ 2 - \Pok\ 



1. 



(4.23) 



This relation together with Eq. (|4.4|1 implies that the power spectrum cannot significantly deviate from C 2 (H2)H\j2tp z 
in the region where /?oo and /3o« are negligibly small. We can also demonstrate the cancellation by explicit calculation 
in the low and high energy limits. For p\r\o\ 1, we have a^ K ss Uo K . Then from Eqs. (|2.20() . 12.21|l . and l|2.32[) . we 
obtain 



dn\a 0ti \ 




2Ktanh(7TK)<iK . , 0/ „ .« / SH 

„ (^ + 1/4)(k» + 9/4) X ^^ U 

nj (k 2 + 1/4) 2 (k 2 + 9/4) x ™ [h^J 



(m « i), 



> i), 



(4.24) 



which gives (pi]o) 2 (£Hi) 2 (5H/Hi) 2 in the low energy regime and (3/ '4)(pr]Q) 2 (5H j Hi) 2 in the high energy regime. 

Comparing these with Eq. (|4.21(l . we see that Vf£{p) cancels P$(p)^. 

To summarize, we have observed that the agreement between the rescaled spectrum P ies (p) and the five-dimensional 
spectrum P§T)(p) is exact up to first order in SHjH\. The agreement is not exact at second order, but we found that 
the correction is not enhanced at any wavelength irrespective of the value of £H\. Just for illustrative purpose we 
show the results of numerical calculations in Fig. [3| and Fig. 0] 
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FIG. 3: Five-dimensional power spectrum of gravitational waves Vzoip) and the four-dimensional one 7 3 4d(p) at low energies 
(Mi = 1CT 2 ) with SH/H! = 10" 3 . These two agree with each other. In this case initial KK modes give negligible contribution. 




FIG. 4: Power spectra of gravitational waves at high energies (£Hi — 10 3 ) with SH/Hi — 1CP 3 . Five-dimensional spectrum 
T , 5^°(p)+T > ^u(p) an d the rescaled four-dimensional one Vr OS {p) agree well with each other (solid lines), while only the zero-mode 
contribution V^u°(p) gives the reduced fluctuation amplitude inside the horizon (dotted line). 

V. DISCUSSION 

In this paper we have investigated the generation of primordial gravitational waves and its power spectrum in the 
inflationary braneworld model, focusing on the effects of the variation of the Hubble parameter during inflation. For 
this purpose, we considered a model in which the Hubble parameter changes discontinuously. 

In the case of de Sitter inflation with constant Hubble parameter H, the spectrum is known to be given by Eq. 1)4. 1|) 
[l3|. It agrees with the standard four-dimensional one (Eq. 1)4. 2f> ) at low energies £H <C 1, but at high energies 
£H ^ 1 it significantly deviates from Eq. I|4.2[) due to the amplification effect of the zero-mode normalization factor 
C(H). One can say, however, the five-dimensional spectrum is obtained from the four-dimensional one by the map 
H^HC(H). 

In a model with variable Hubble parameter, gravitational wave perturbations are expected to be generated not only 
from the "in-vacuum" of the zero mode but also from that of the Kaluza-Klein modes. Hence, it is not clear whether 
there is a simple relation between the five-dimensional spectrum and the four-dimensional counterpart. Analyzing 
the model with a discontinuous jump in the Hubble parameter, we have shown that this is indeed approximately the 
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case. More precisely, if the squared amplitude of four-dimensional fluctuations is given by (/i/2-7rMpi) 2 , we transform 
(h/2irMpi) 2 into C 2 (h)(h/2irMp{) 2 , then the resulting rescaled spectrum exactly agrees with the five-dimensional 
spectrum Vsoip) up to first order in 8H/H1. At second order 0(SH / H) 2 the agreement is not exact, but the 
difference is not enhanced at any wavelength irrespective of the value of IH\. Hence, in total, the agreement is not 
significantly disturbed by the mismatch at second order. As a non-trivial point, we also found that the initial KK 
mode vacuum fluctuations can give non-negligible contribution to the final zero- mode states at second order. 

One may expect that the power spectrum of gravitational waves in the braneworld model would reflect the charac- 
teristic length scale corresponding to the curvature (or "compactification" ) scale of the extra dimension I. However, 
our analysis showed that the resultant power spectrum does not depend on the ratio of the wavelength of gravitational 
waves to the bulk curvature scale, p\ij \£H. 

Here we should mention the result of Ref. [l8| that is summarized in Appendix ^ Their setup is the most violent 
version of the transition Hi — > Hi — 0. If the wavelength of the gravitational waves is much longer than both the 
Hubble scale and the bulk curvature scale, the power spectrum is given by Eq. I|A4(1 and obviously it is obtained from 
the four-dimensional counterpart by the map H 1— ► HC(H). However, if the wavelength is longer than the Hubble 
scale but much smaller than the bulk curvature radius, the amplitude is highly damped as is seen from Eq. (| A7|) and 
the map H t— > H C (H) does not work at all. This damping of the amplitude can be understood in the following way. 
In the high energy regime IH 3> 1, the motion of the brane with respect to the static bulk is ultrarelativistic. At the 
moment of transition to the Minkowski phase, the brane abruptly stops. Zero-mode gravitons with wavelength smaller 
than the bulk curvature scale can be interpreted as "particles" traveling in the five dimensions. These gravitons make 
a "hard hit" with the brane at the moment of this transition, and get large momenta in the fifth direction relative to 
the static brane. As a result, these gravitons escape into the bulk as KK gravitons, and thus the amplitude (|A7() is 
damped. On the other hand, in our model the change of the Hubble parameter is assumed to be small, and hence such 
violent emission of KK gravitons does not happen. If this interpretation is correct, the mapping rule h 1— > hC(h) will 
generally give a rather good estimate for the prediction of inflationary braneworld models as far as time variation of 
the Hubble parameter is smooth. If we can confirm the validity of this prescription in more general cases, the analysis 
of gravitational wave perturbations will be simplified a lot. We would like to return to this issue of generalization in 
a future publication. 
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APPENDIX A: PARTICLE CREATION WHEN CONNECTED TO MINKOWSKI BRANE 

Here for comparison with our results we briefly summarize the results obtained by Gorbunov et al. [T§ | focusing 
on the power spectrum of gravitational waves. Their method is basically the same that we have already explained in 
the main text. They considered the situation that de Sitter inflation on the brane with constant Hubble parameter 
H suddenly terminates at a conformal time 77 = 770 , and is followed by a Minkowski phase. The power spectrum of 
gravitational waves is expressed in terms of the Bogoliubov coefficients as 

V5d(p) = W (|Q (w°) 2 ( x + 2 \M 2 + 2 J ^IA)k| 2 ) , (Ai) 

where |/3oo| 2 arl d |/3o K | 2 are the number of zero- mode gravitons created from initial zero- mode and KK modes, respec- 
tively. At super Hubble scale (p\t]q\ <C 1) we can neglect the first term in the parentheses, which corresponds to the 
vacuum fluctuations in Minkowski space. 

According to when p\r] \£H <C 1 (i.e., the wavelength of gravitational wave (p/a)^ 1 is much larger than the 
bulk curvature scale i) and p\rjo \ <C 1, the coefficients are given by 

l/^l 2 * ggl (A2) 
dKlM ~ { P \vo\tH (IH»1). (A3) 
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One can see that the contribution from initial KK modes is suppressed irrespective of the expansion rate £H . Therefore 
the power spectrum is evaluated as 

Equation (|A4|I is half of the power spectrum on the de Sitter brane (Eq. (|4.1|l ). and this result can be understood as 
follows. The amplitude of fluctuations at the super Hubble scale stays constant. After the sudden transition from the 
de Sitter phase to the Minkowski phase, the Hubble scale becomes infinite. Therefore, those fluctuation modes are 
now inside the Hubble horizon, and they begin to oscillate. As a result, the mean-square vacuum fluctuation becomes 
half of the initial value. 

On the other hand, when p\rjo\£H ^s> 1 and p\rjo\ <C 1 (these conditions require IH ^> 1), the Bogoliubov coefficients 
are given by 

|/300| 2 « C2(H L ]rr^ (A5) 



ipvo) 2 ipnoiH) 

As before, the contribution from initial KK modes is negligible, and that from | /?oo 1 2 dominates the power spectrum, 

C 2 (H) ( H\ 2 4 

P5D "^rUJ j^my- (A7) 

One can see that the spectrum is suppressed by the factor 4/ (pi]o£H) 2 . 

APPENDIX B: CALCULATIONS OF THE INNER PRODUCT 

We derive Eqs. I|2.31[) and l|2.32[l by calculating the inner product (|2.3U|) . Expanding Eq. (|2.25|) in terms of ch, we 
have 

fj = r ? + eff?7ocosh£-4^(sinh£) 2 /(2r7) + --- , (Bl) 
I = £ - £ff % sinh £_/t] + e 2 H ^ cosh £ sinh ^/rj 2 + ■■■ , (B2) 

which reduce to the following forms at the infinite past (rj — > -co), 

fj-r) = e//?7ocosh^, (B3) 
t-t = 0. (B4) 

Then, Uqo is evaluated at 77 = — 00 as 



U 00 = -2ieQ ^ 2(s f h0 3 {fc ] dA +) ~ ti +) d n &- ) 

jj rca e -ie H pr) cosh J 



~ M^l\ e ^HPn co S ^ B . 2 £ 2 H 2 C 2 (H 1 ) x 



roo 

xe i f -Hpr)o coshes / ^ 



1 - ie H pr)oCOsh£, - i(e H Wo) 2 (cosh^) 2 



l (B s (sinh0 3 

where we expanded the integrand with respect to en in the last line. Integrating each term, we finally obtain 



TT HlCjH-l) —je H pr) coshfs 

00 ~ HxG{H x ) 



1 - ie H pi] (C 2 (iJi) - coshes) - i(e if p?7o) 2 (sinh^ B ) 2 



(B5) 



(B6) 



Note that the condition e#p|?7o| coshes (w p\rjo\SH/H) < 1 is required in order to justify the expansion of the 
exponent. Because the integral is saturated at £ sa we do not have to worry about the validity of the expansion 
for large cosh£. 
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The explicit form of Uo K is needed up to the first order in eh- A similar calculation leads to 



U 0K 



-2U 

' Cb 77 2 (sinh^) 
= -2i£ 3 ■ r 1/2 C(H 2 ) 



(-) 
o 



-3/2 roo 



2p y/2p J iB ?7 2 (sinhO : 



r\ — ► — oo 

-MO 



fj--\ e - ipfi • (-1 - ipr 1 )e im+i7r ' i - (-n) e ipj?+i7r / 4 • (-ipr])e'^ 
P) ' . 

roc p -i£HPVo cosh^ B 

-2^HH 2 CiH 2) j g (8inh03 Xw (fl. 



(B7) 



The integral in the last line, which we call /, can be calculated as follows. Again, expanding the integrand in terms 
of 6h, we have 



I : 



(sinhO 3 



(1 - ie H pm cosh f) 



(B8) 



Let us consider the first term in the parentheses. The spatial wave function \k satisfies (sinh^)~ 3 x K = — (k 2 + 
9/4)~ 1 d^[(a.nh^)~ 3 d^XK] (Eq. J^TOJl) with the boundary condition O^Xk^b) = 0- Therefore, together with the 
behavior at infinity, (sinh£) _3 (?£X K ~ (sinh^) _3 9^(sinh^) 3 / 2 — > 0, we find that the integral of the first term vanishes. 
Then, using the integration by parts twice, we have 



k + — ] / w itHPVo / dt; cosh £ 



d_ 



d 



= -lenprio 



1 



(sinh£) 3 <9£ 



itHPVo 
itHPVo 



(sinh^ B ) 



(sinh£) 2 <9£ 

p -2coshe _ 



2/, 



(B9) 



from which we can evaluate J7o«- Note that the approximation is valid when p\r)o\5H/Hi -C 1. 
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